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We show that an appropriate choice of the potential parameters in one-dimensional quantum 
systems allows for unity transmission of the tunneling particle at all incident tunneling energies, 
except at controllable exceedingly small incident energies. The corresponding dwell time and the 
transmission amplitude are indistinguishable from those of a free particle in the unity-transmission 
regime. This implies the possibility of designing quantum systems that are invisible to tunneling by 
a passing wave packet. 
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I. INTRODUCTION 

The design and construction of one-dimensional artifi- 
cial quantum structures at nanometric scales has opened 
a new realm of possibilities on the investigation of fun- 
damental properties of quantum mechanics [T]. One of 
these properties is tunneling which represents one of the 
paradigms of quantum mechanics. As discussed in quan- 
tum mechanics textbooks, tunneling of a particle of a 
given energy through a potential barrier yields in general 
partial transmission. Full transmission is exhibited in 
resonant tunneling systems, which at least are formed by 
two barriers with a well in between. There, unity trans- 
mission may be achieved at some specific energies, the 
so called resonance energies [2]. This yields, however, a 
time delay with respect to free propagation, that is pro- 
portional to the inverse of the resonance energy width 
OS], and hence it allows one to distinguish the tunnel- 
ing particle from one evolving freely. The issue of total 
transparency of a tunneling particle by a potential along 
the full energy range has attracted attention over the 
years. It has been addressed within different frameworks: 
inverse scattering theory [5], supersymmetric quantum 
mechanics |6j, Darbox transformation approach [7] and 
group-theoretical approaches 0. These works refer to 
a number of exactly solvable potentials, usually named 
reflectionless or transparent potentials, for which the re- 
flection amplitude vanishes identically, while the trans- 
mission amplitude has modulus 1 for all incident energies 
E including the threshold energy value E — 0. A well- 
known example is the Poschl- Teller (P-T) potential well, 
which for very specific values of the potential parameters 
attains unity transmission at all energies |^. However, 
transparent potentials have escaped, to the best of our 
knowledge, experimental verification and are mainly of 
interest in mathematically oriented studies. A possible 
reason is that transparency in these potentials is tightly 
bound to the functional dependence of the potential. 

Here we investigate to what extent one may design po- 
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tential profiles in one dimension (ID) that, in addition to 
being totally transparent to a tunneling particle, cannot 
be detected by interference experiments. Our motivation 
is purely theoretical and would lead to the possibility of 
designing invisible quantum systems. Our approach rests 
on analytical properties of the outgoing Green's function 
of the system in the complex momentum plane that hold 
provided the potential vanishes beyond a distance and 
the transmission is a coherent, elastic process. These 
analytical properties consist of having a bound or an an- 
tibound pole very close to the energy threshold and all 
other poles far away and overlapping among themselves. 

We find that one may design potential profiles in ID 
that possess two regimes for transmission of incident 
monochromatic energy particles. In one regime, occur- 
ring at very small controllable energies close to the energy 
threshold, i.e., a very small fraction of the potential bar- 
rier height, the transmission coefficient rises sharply from 
zero to unity, and in the other regime, which involves 
the rest of tunneling energies and energies extending up 
to several times the potential barrier height, the parti- 
cle attains essentially unity transmission. We show that 
in the unity-transmission regime the transmission phase 
has a vanishing value, which implies indeed that interfer- 
ence experiments cannot detect the scattering potential 
and, in addition, that the dwell time, which provides the 
relevant time scale for the tunneling process, is indistin- 
guishable from that of a free particle. As a consequence 
of the above considerations we find that in the unity- 
transmission regime these systems are indeed invisible to 
a tunneling particle. Moreover, since in time domain, 
extremely small energies correspond to very long times, 
we obtain that these systems are essentially invisible to 
tunneling by an incident pulse or wave packet. We shall 
refer to these systems as invisible systems. 

It is worth noticing that here invisibility refers to a 
different process from studies that involve the design of a 
cloak surrounding a system that then becomes invisible 
to light |10j or with approaches in the quantum domain, 
refer to as quantum cloaking, where a system is sur- 
rounded by a cloak to become invisible to matter waves 
at certain incident energies in two dimensions (2D) and 
three dimensions (3D) [TTJ [T^]. These approaches are 
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based on ideas from transformation optics and refer only 
to a time-independent description. We do not surround 
a system with a cloak but rather we design systems that 
become invisible to matter waves and consider both the 
energy and time domains. 

Contrary to transparent potentials where full trans- 
mission is tightly bound to the functional dependence 
of these potentials, the potentials considered here are ro- 
bust against some variation in the functional dependence 
of the potential profile. One may consider rectangular 
or continuous shapes formed by distinct combinations of 
barriers and wells. 

This work is organized as follows. In Sec. II we con- 
sider the resonance formalism and the relationship be- 
tween the transmission amplitude and the distribution of 
its complex poles. Section III deals with invisible systems 
through a number of subsections that discuss, respec- 
tively, the rectangular barrier and the Poschl- Teller po- 
tential, multibarrier systems, the dwell time, and wave- 
packet scattering. Finally, section IV, gives the conclud- 
ing remarks. 

II. TRANSMISSION AMPLITUDE AND 
COMPLEX POLES 

Let us consider a particle of mass m and energy E 
impinging, from x < 0, on a quantum structure char- 
acterized by a potential profile V{x) of length L, i.e., 
V{x) = outside the region < x < L. As is 
well known, the solution to the Schrodinger equation of 
the problem may be written for x < as ^l),^{x,t) = 
exp(ikx) + r(k) exp(—ikx), and for x > L as ^-^{Xjt) = 
t(fc) exp(ifca:), where T{k) and t(fc) stand, respectively, 
for the reflection and transmission amplitudes. It is 
convenient to write the transmission amplitude t(fc) in 
terms of the outgoing Green's function of the problem, 
G'^{x,x';k) [13^, namely, 

t(fc) = 2iA:G+(0,L;fc)e-''=-^, (1) 

where k = [2mE]^^^/n. The reason is that this allows one 
to obtain a representation for the transmission amplitude 
as an expansion involving the poles and residues of the 
outgoing Green's function to the problem. This proce- 
dure is in fact numerically equivalent to standard nu- 
merical calculations such as the transfer-matrix method 
[T]. However, it yields a deeper physical insight by es- 
tablishing a link between the analytical properties of the 
outgoing Green's function on the complex k plane and 
the behavior with energy of the transmission phase and 
the transmission coefficient. 

It is well known that the function G'^{x,x';k), and 
hence the transmission amplitude t(fc), possesses an infi- 
nite number of complex poles fc„, in general simple, dis- 
tributed on the complex k plane in a well-known man- 
ner (31 1^ . Purely positive and negative imaginary poles 
kji = ijn correspond, respectively, to bound and anti- 
bound (virtual) states, whereas complex poles are dis- 
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FIG. 1: (Color online) (a) Transmission coefficient as a func- 
tion of tfie energy in units of tlie potential fieight Vq for 
two double-barrier systems. The exact numerical calculation 
(solid line) is reproduced exactly by Eq. Q using A*' = 500 
poles. Both systems have the same barrier height, Vb = 0.2 
eV, and no well depth, Uo = 0; their well widths are twice the 
barrier widths, i.e., w — 26. One of them (dotted line) has 
b = 4.0 nm and exhibits two well-defined resonances along 
the tunneling region, whereas the other (dashed line), b = 0.4 
nm, shows no resonance structure at all. (b) Distribution of 
several complex poles of the transmission amplitude in the 
l3 = kL plane for the potentials in (a): for b — 4.0 nm (cir- 
cles) and for b = 0.4 nm (stars). One sees that by diminishing 
the values of b, and hence of w, one goes from a system with 
sharp resonances (poles very close to the real axis) to one with 
no resonances at all (all poles far away from the real axis). 
See text. 



tributed along the lower half of the k plane. They may be 
calculated by using iterative techniques as the Newton- 
Raphson method [H]. The outgoing Green's function 
G^(0, L; k) may be expanded as an infinite sum in terms 
of its poles [131 IE] • We have found recently that the ex- 
pansion of G~^{0, L; k) exp{—ikL) has better convergence 
properties. It yields the expansion for the transmission 
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amplitude 



t{k) = 2ik 



iknL 



n— — oc 



(2) 



where r„ follows from the residue of G"'"(a;,a;'; k) at the 
pole kn [ij 113] . The position of the poles fc„ on the com- 
plex k plane is a function of both the parameters of the 
potential and the mass of the particle. Consequently, 
by varying these parameters the poles follow trajectories 
along the k plane. 

For a given combination of rectangular barriers and 
wells, we denote, respectively, the barrier heights and 
depths by Vq and — C/q, measured in eV, and the rectan- 
gular barrier and well widths by h and w, measured in 
nm. This is sufficient to characterize a variety of possible 
combinations of rectangular barriers and wells, as the 
barrier-well (BW), the barrier-well-barrier (BWB), the 
well-barrier-well (WBW) systems, and so on. In order 
to exemplify the above considerations and the relation- 
ship of pole distributions with the behavior of the trans- 
mission coefficient as a function of energy, which follows 
from Eq. we consider two double-barrier tunneling 
systems (BWB) with parameters typical of semiconduc- 
tor tunneling structures pj, as indicated in Fig. [l] In all 
calculations the effective electron mass is taken as that of 
GaAs, i.e., m = 0.067 me, with rrie as the free-electron 
mass. Figure [ija) provides a plot of the transmission 
coefficient as a function of energy in units of the poten- 
tial height, which is the same for both systems. In both 
systems the well depths are zero and the well widths are 
twice the barrier widths. In system 1 (dotted line), the 
barrier and well widths are ten times larger than in sys- 
tem 2 (dashed line). One sees that system 1 exhibits 
two well-defined resonances along the tunneling region, 
whereas system 2 exhibits no resonances at all. The 
above behavior of the transmission coefficient reflects it- 
self in the distribution of the complex poles A:„ = /z„ — if„ 
of the corresponding transmission amplitude, shown in 
Fig. [l]^b). In the case of system 1 (circles), there appear 
two complex poles very close to the real f3 = kL axis and 
one may follow well known arguments to show that each 
of them yields a Lorentzian or Breit-Wigner analytical 
expression for the transmission coefficient near resonance 
energy [16l [17] . From the third pole onward the width of 
the poles increases steadily and one sees that the trans- 
mission coefficient eventually approaches unity. In the 
case of system 2 (stars), the poles are all situated away 
from the real axis, except for an antibound pole situated 
no far from the threshold value. We shall see below how 
important poles near threshold are for invisibility. Notice 
that since complex poles obey, from time-reversal invari- 
ance considerations |16| . the relationship fc_„ — — fc*, 
only poles seated on the fourth quadrant of the /3 plane 
have been depicted. 




FIG. 2: (Color online) Potential profiles of a two-double- 
barrier rectangular (2BWB) system (dashed line) and a two- 
double-barrier P-T potential (dotted line). 



III. INVISIBLE SYSTEMS 

A. Rectangular barrier and Poschl- Teller potentials 

Recently it has been shown that total transparency 
of a very thin single-barrier rectangular potential at all 
except very small energies follows from a distribution of 
poles that consists of an antibound pole seated very close 
to fc = and all other complex poles away from the real k 
axis and overlapping with each other [T5]. The antibound 
pole ka may written as [T5]. 
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A similar situation holds for the P-T barrier potential 
W{x) = Vol cosh^{x/d) [9]. Here, we may also denote, 
respectively, the corresponding barrier height or depth 
by Vq or —Uq and the barrier or well widths by the pa- 
rameters c?b or d^,. The transmission amplitude for the 
P-T barrier potential reads 



t(fc) 



sinh(7rfcd)e 



sinh(7rfc(i) -f zcos[(7r/2)-yi — ry] ' 



(4) 



with <j) as & phase and -q = SmUod'^/h^. Equation ^ has 
poles at the zeros of its denominator. For 77 ^ 1, the 
P-T potential has an antibound pole very close to fc = 0, 

namely, at 
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[2mVo]d 
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which resembles that for the thin rectangular barrier po- 
tential written above. Clearly for a P-T well, where the 
potential parameter is negative, i.e., — C/q, a similar re- 
lationship holds for a bound state 7^. The above an- 
alytical behavior is different from the well known total 
transparency of a single P-T well at all energies, includ- 
ing E = 0, which occurs for {1 + rj) = (2n -1-1)^, with 
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FIG. 3: (Color online) Distribution of the complex poles of the 
transmission amplitude in the /3 = kL plane for several poten- 
tial profiles: BWB (triangles), 2BWB (stars), 5BWB (dots), 
and quadruple-barrier Poschl- Teller potential (circles). Each 
potential is characterized by having a bound or an antibound 
pole very close to the threshold P — (see inset) and all other 
poles overlapping and away from the real /3 axis. See text. 



n — 0,1, 2, ... [9,. In this case 77 may be quite large and 
the corresponding outgoing Green's function has, as only 
singularity, a pole at fc = 0. The above results for near 
energy threshold bound or antibound poles suggest to 
look for a similar behavior in systems formed by different 
combinations of barriers and wells for either rectangular 
or P-T potential profiles. In the case of P-T potentials 
this necessarily introduces a cutoff in the potential tails 
and hence the analytical properties of the transmission 
amplitude become analogous to that of rectangular po- 
tentials. 



B. Multibarrier systems 

For both rectangular and continuous potential pro- 
files one may consider different combinations of BWB or 
WBW systems to form, for example, chains of these sys- 
tems, as the quadruple-barrier system (2BWB) formed 
by two BWB systems separated by a distance h, etc. Fig- 
ure [2] illustrates the potential profiles for a two-double- 
barrier rectangular potential and a two-double-barrier P- 
T potential. 

Figure [3] provides examples of these pole distribu- 
tions on the P = kL plane, for a number of systems: 
a BWB (triangles), a quadruple-barrier 2BWB (stars), 
a ten-barrier 5BWB (dots), and a quadruple-barrier P- 
T potential (circles). For all the rectangular potential 
systems, we consider also parameters typical of semi- 
conductor heterostructures: b = 0.4 nm, w = 0.8 nm, 
h = 0.8 nm, and also Vq = |C/ol = 0.12 eV, except 
for the 5BWB system, where the depth of the second 



FIG. 4: (Color online) The transmission coefficient T{E) as a 
function of energy in units of the potential height Vb is calcu- 
lated, using Eq. Q, for some of the systems considered in Fig. 
|3] 2BWB (dotted) and 5BWB (short-dashed line) for rectan- 
gular barrier-well potentials, and quadruple-barrier P-T po- 
tential (dashed line) . Also shown is an exact calculation for a 
quadruple-barrier system (2BSB) with the same parameters 
as that of 2BWB except that the well depths U = G (short 
dots). In this case the transmission is not unity along the 
tunneling region. All calculations are reproduced exactly by 
rmmerical calculations of T{E), as exemplified for the 5BWB 
system (solid line). See text. 



and fourth wells is U = -0.113 eV. For the P-T po- 
tential we choose = |C^o| = 0.12 eV, 4 = 0.0709 
nm, and d.^, = 0.1399 nm. The effective electron mass is 
taken also as in the examples considered in Fig. [T] i.e., 
m/nie = 0.067. It is worth noticing that in all examples 
i^n > T^/L, which establishes a scale for the distance from 
the real k axis of the overlapping complex poles, which 
fulfill (^n-i-i — Mn) n/L. The inset shows a zoom of the 
positions of bound and antibound poles close to fc = 
for the above systems. The values of bound or antibound 
poles may be controllable by choosing appropriately the 
parameters of the potential, as the 5BWB potential ex- 
emplifies. Notice that in order to obtain values for the 
bound or antibound poles so close to the threshold, avoid- 
ing extremely small values of the barrier widths, it seems 
necessary that the well depths Uq differ from zero. 

From an analytical point of view the above results for 
the distribution of the complex poles suggest that the 
outgoing Green's function in these systems is governed, 
similarly to the transparent rectangular barrier |18j . by 
the purely imaginary pole seated close to the threshold 
fc 0, namely. 



G+iO,L;k) 



1 



ikL 



2i{k - i^q) 



(6) 



where q = a or b refers, respectively, to antibound or 
bound pole and l/2i follows from the residue at the 
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imaginary pole kq = ijq |18j . We have verified numeri- 
cally the validity of the above value of for the distinct 
systems considered. Substitution of the expression for 
G+(0, L; k) given by Eq. ^ into Eq. Q yields 



t(fc) 



1 - i^q/k' 



(7) 



where we have used exp{—ikqL) « 1. Notice that /c^ « 
implies that the modulus of t is very close to unity and 
that its corresponding phase 9 « Jq/k is close to zero 
except at very small values of k and hence of energy. 
It is worth noticing that the expression for G^{0,L;k), 
given by Eq. (|6|, exhibits a singularity very close to 
k = 0, which resembles the singularity at fc = of the 
free outgoing Green's function, 



G+(0,i;fc) 



1 



(8) 



It follows from Eq. ([t]) that the transmission coefficient 
reads. 



T{E) = \t{E)\'' 



1 



l+Eq/E' 



(9) 



where Eq ~ 



{n^/2m)-fl. Figure [i] exhibits a plot of 
T{E) as a function of energy in units of the poten- 
tial height Vo, for several of the systems considered in 
Fig. |3] 2BWB (dotted) and 5BWB (short-dashed line), 
for rectangular barrier-well potentials and a quadruple- 
barrier P-T potential (dashed line) as that depicted in 
Fig. [2] The corresponding values of Eq for the these 
potentials are, respectively, -E2BWB = 8.68 x 10~^ eV, 
i^sBWB = 1-67 X 10-^ eV, and E^p-t = 6.19 x 10^1° 
eV. It might be of interest to compare the above values 
of Eq, for multibarrier rectangular systems, with that of 
a single rectangular barrier. This follows by substitu- 
tion of Eq. ([5]) into the above expression for Eq to give 
Eq — [{2m/h )Vq /A\L'^. For example, for a barrier of 
both, with similar height (i.e., Vq — 0.12 eV) and effec- 
tive mass (i.e., m/me — 0.067), a value of Eq ^ 10^^ 
eV would require a width h = L = 0.012 nm and for 
Eq ~ 10~® eV, 5 = L = 0.0012 nm. The above values for 
L are extremely small. The widths of barriers and wells in 
multibarrier systems along the unity-transmission regime 
are much larger than for a single-barrier system. A sim- 
ilar situation holds regarding P-T potentials. 

The calculations using Eq. ^ are indistinguishable 
from the corresponding exact numerical calculations us- 
ing the transfer-matrix method 1 . The differences 
among the distinct systems are only appreciable at very 
small energies. Figure |4] exhibits also the transmission 
coefficient for a quadruple-barrier potential with poten- 
tial depths U—Q, 2BSB, (short dots). This case is similar 
to the BSB system (dashed line) presented in Fig. [l] Al- 
though it possesses both overlapping complex poles and 
an antibound pole, the energy of this antibound pole, 
^2BSB = 4.18 X 10~^ eV, is not sufficiently close to the 
energy threshold to exhibit unity transmission along the 
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FIG. 5: (a) Transmission contour as a function of the energy 
E and the parameter V = Vo = \U\, for quadruple-barrier 
systems 2BWB, and all other parameters fixed, for V > Q, 
and quadruple-well systems 2WBW, for V < Q. (b) Trans- 
mission contour as a function of the energy E in units of the 
potential height Vb for a quadruple-barrier system 2BWB vs 
log]^g(m/me), where m is an effective mass and stands for 
the free-electron mass. See text. 



tunneling region. Notice that it is several orders of mag- 
nitude larger than the values for the other systems. 



1. Robustness 

The phenomenon of invisibility is robust against some 
variation in the values of the potential parameters of 
the invisible system. Within certain limits the variation 
in effective masses, barrier heights, well depths, barrier 
widths, and well widths, either for rectangular or con- 
tinuous shapes as the P-T potentials, keeps the system 
invisible. To exemplify this. Fig. [sj^a) exhibits a con- 
tour plot for the transmission coefficient as a function 
of the energy E and the parameter V = Vq = \U\ ior 
rectangular quadruple-barrier systems 2BWB, where all 
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the other potential parameters have the same values as 
given previously. Notice that for negative values of V the 
above systems become quadruple-well systems 2WBW. 
It is also worth noticing that along the "invisibility win- 
dow", the 2BWB and the 2WBW systems are indistin- 
guishable from each other. The plot for the contour of 
the transmission coefficient T{E) considers the range of 
values 0.5 < T{E) < 1, which is the range employed for 
resonance transmission. There is a range of values of V 
around V = 0, the free case, that correspond to full trans- 
parent systems. One may also consider a similar varia- 
tion regarding the barrier and well widths, and again 
within certain limits, full transparency remains robust. 
Figurejsjb) exhibits the effect of the variation in the effec- 
tive mass m for the rectangular quadruple-barrier system 
2BWB discussed above. This figure displays a contour 
plot for the transmission coefficient, where logiQ^m/me) 
varies in a broad range of values for different incidence 
energies in units of the barrier height Vq. The value 
used in the previous calculations, m/nie — 0.067, which 
corresponds to a GaAs quadruple-barrier P-T potential 
(circles), yields log]^Q(0.067) = —1.1739, which clearly 
falls within the invisibility regime. The same occurs for 
m/nie = 0.1, where log]^Q(0.1) = —1.0, a value commonly 
used for GaAsAl barriers. Notice that as the effective 
mass increases, the system eventually ceases to be invis- 
ible and may exhibit a resonance structure. 



C. Dwell time 

Let us now investigate the dwell time in these sys- 
tems. The dwell time is defined as 



1 

Jo 



mx,E)\'^dx, 



(10) 



where Jq = hk/m stands for the incoming flux. This 
quantity measures the amount of time that the incident 
particle spends within the internal region. One may write 
it in units of tq = L/ Jq, the time it takes to a free particle 
to traverse the distance L, and express it as [20| |2T]. 



Td 



\i:{x,E)fdx = T+]r Te + R(f> 

±j 



kL 



■ sm ( 



where R stands for the reflection coefficient, 9 and (j> refer, 
respectively, to the so called transmission and reflection 
times, the dot representing the derivative with respect to 
k of the phases 9 and 4> of the corresponding transmission 
and reflection amplitudes t(fc) and r(fc). Figure [6] yields 
a plot of Td{E) in units of tq vs E for several systems: 
2BWB (solid hue) and 5BWB (dotted line ), with pa- 
rameters as considered above, and lOBWB (dashed line), 
which is formed by two 5BWB systems separated also 
by a distance h = 0.8 nm. The lOBWB system has a 
length of L = 23.2 nm and possesses an antibound pole at 
— 1.09118 X 10"'^ nm~^. The exact numerical calculation 




FIG. 6; (Color online) The integral expression for the dwell 
time Td in units of tq = L/ Jq, as a function of E/Vo, is eval- 
uated numerically for the systems 2BWB (solid line), 5BWB 
(dotted line), and lOBWB (dashed line). Notice that except 
at very small energies Td is very close to tq. 



of Td{E) is obtained by integrating the probability den- 
sity along the internal region of the potential using the 
transfer-matrix method. One sees that in all cases, ex- 
cept at very small energies, Td{E) is very close to tq{E). 
The above result implies that the sum of the last two 
terms on the right-hand side of Eq. ( 11 ) adds to a van- 



ishing contribution, although each term by itself may not 
be small. One may conclude that in these systems, ex- 
cept at very small energies, the time that the tunneling 
particle spends along the internal region of the potential 
is indistinguishable from that of a free evolving particle. 
These systems might be used to make a comparison of the 
different definitions for tunneling times, which remains a 
long-debated and controversial subject 



D. Wave packet scattering 

The above discussion refers to monochromatic waves. 
Let us now consider the tunneling of a Gaussian wave 
packet on these systems. The initial wave packet ip{x, 0) 
is represented by 



V'(a;,0) = Ae-(^- 



(12) 



satisfying the condition |a::o|/2(T > 1, which guarantees 
that the tail of the Gaussian wave packet is very small 
near the interaction region < a; < L, and may be 
solved analytically [24^. We could make a comparison 
along the transmitted region of a wave packet evolv- 
ing freely, \il}f{x,t)\^, with the wave packet that tun- 
nels through the system, |^(a;,t)p. However, even if 
\'ipf{x,t)\^ = \^{x,t)\'^ as a function of time, it is not 
sufficient to conclude that the system is invisible, since 
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Re{tl^*j:{x,t)il:{x,t)} = cos{0)\^|:*{x,t)\^. Thus if cos(6') = 
1, we may conclude that the system is invisible to the 
tunneling wave packet beyond any doubt. Figure [7] il- 
lustrates that this is indeed the case. It yields a com- 
parison of ^{x,t) — Re t)} (solid hne) for 
the system 2BWB discussed above with the correspond- 
ing free time evolution wave packet pf(x,t) = \'ipf{x,t)\'^ 
(dots) as a function of time in units of to = {x — xo)/vo. 
One sees that both solutions are indistinguishable from 
each other. The inset exhibits a similar comparison for 
the system 2BSB, whose only difference with the system 
2BWB is that the well depths are zero, and hence it does 
no exhibit unity transmission along the tunneling region 
as shown in Fig. |6] 



FIG. 7: (Color online) Comparison of two initially iden- 
tical Gaussian wavepackets, one evolving freely (dots), de- 
scribed by pf{x,t) = \'ipf{x,t)\'^, and the other tunneling 
through a 2BWB system (solid line). Here we calculate 
^{x, t) — Re {ip*f{x, t)ip{x, t)} = cos{0)\tlj}{x, t)\^ to show that 
there is no phase dependence on the transmitted phase 6 and 
hence the system is invisible. In contrast, the inset shows that 
a similar calculation for the 2BSB system considered in Fig|4] 
exhibits a phase dependence. The initial parameters of the 
Gaussian are: a = 0.5 nm and xo — —5.0 nm with an energy 
E = 0.06 eV (half the barrier height). The above quantities 
are calculated at x — 100.0 nm as a function of time in units 
of to = {x - xo)/vo. 



there might exist a dependence on the transmitted phase. 
Indeed if we write 

,P{x,t) = e'',Pf{x,t), (13) 

in order to exhibit a possible phase dependence 
it is more convenient to compare \'ipf{x,t)\'^ with 



IV. CONCLUDING REMARKS 

In summary, we predict the possibility of designing ar- 
tificial quantum systems in ID that are invisible to a 
passing wave packet. Hence the system becomes unde- 
tectable by matter waves. Although our examples refer 
to rectangular and P-T multibarrier systems, they are 
of a general nature in quantum physics and may also be 
considered in other artificial systems as ultracold atoms 
in optical lattices. Our results depend on general analyt- 
ical properties of the transmission amplitude for coherent 
processes and may open the way to the design, experi- 
mental scrutiny, and applications of these quantum sys- 
tems. 
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